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We systematically study topological phases of insulators and superconductors (or superfluids) in 
three spatial dimensions (3D). We find that there exist 3D topologically non-trivial insulators or 
superconductors in five out of ten symmetry classes introduced in seminal work by Altland and 
Zirnbauer within the context of random matrix theory, more than a decade ago. One of these is 
the recently introduced Z2 topological insulator in the symplectic (or spin-orbit) symmetry class. 
We show there exist precisely four more topological insulators. For these systems, all of which 
are time-reversal invariant in 3D, the space of insulating ground states satisfying certain discrete 
symmetry properties is partitioned into topological sectors that are separated by quantum phase 
transitions. Three of the above five topologically non-trivial phases can be realized as time-reversal 
invariant superconductors, and in these the different topological sectors are characterized by an 
integer winding number defined in momentum space. When such 3D topological insulators are 
terminated by a two-dimensional surface, they support a number (which may be an arbitrary non- 
vanishing even number for singlet pairing) of Dirac fermion (Majorana fermion when spin rotation 
symmetry is completely broken) surface modes which remain gapless under arbitrary perturbations 
of the Hamiltonian that preserve the characteristic discrete symmetries, including disorder. In par- 
ticular, these surface modes completely evade Anderson localization from random impurities. These 
topological phases can be thought of as three-dimensional analogues of well known paired topo- 
logical phases in two spatial dimensions such as the spinless chiral (p x ± ip y )-wa,ve superconductor 
(or Moore-Read Pfaffian state). In the corresponding topologically non-trivial (analogous to "weak 
pairing") and topologically trivial (analogous to "strong pairing") 3D phases, the wave functions 
exhibit markedly distinct behavior. When an electromagnetic U(l) gauge field and fluctuations of 
the gap functions are included in the dynamics, the superconducting phases with non-vanishing 
winding number possess non-trivial topological ground state degeneracies. 



I. INTRODUCTION 

Quantum states of matter are characterized not only 
by the structure of the enjaray spectrum but also by the 
nature of wave functions .lilttU Of particular importance 
are topological properties of wave functions, i.e., proper- 
ties that are invariant under small adiabatic deformations 
of the Hamiltonian. A classic example of such topologi- 
cal characteristics is the quantized Hall conductivity o~ xy 
in the integer quantum Hall effect (IQHE), which oc- 
curs at low temperature and high magnetic field in two- 
dimensional (2D) electronic systems with broken time- 
reversal symmetry (TRS). The transverse (Hall) con- 
ductivity a xy , arising from topologically protected edge 
currents, can be intesBreted as an integer Chern num- 
ber (TKNN integer) ,M a quantized topological invari- 
ant which characterizes the different topological ground 
states. 

In the past few years, it has been realized that 
topological phases supporting topologically protected 
states appearing at the sample boundaries can also ex- 
ist in two- and three-dimensional time-reversal invari- 



ant 
field 



lias 



aence of an external magnetic 
]B These topological states occur 
in certain materials with a bulk band gap generated 
by strong spin-orbit interactions and are known as Z2 
topological insulators. Unlike the integer quantum Hall 
states, the two-dimensional version of the Z2 topological 



insulator, which has been dubbed "quantum spin Hall" 
(QSH) state, does not carry any net charge current along 
the edges. Instead, when a U(l) part of the SU(2) spin 
rotation symmetry is conserved, electrons with opposite 
spin propagate in opposite directions, .-athich gives rise 
to a quantized spin Hall conductance.OEJ Remarkably, 
the topological order of the QSH state survives even un- 
der a (small) breaking of the full spin-rotation symmetry. 
Consequently, the QSH insulator cannot be character- 
ized by a quantized spin Hall conductivity. Rather, as 
shown by Kane and Mele in Ref. [?], there is a Z2 topo- 
logical invariant that classifies the topological properties 
of the QSH states in a similar way as the Chern num- 
ber does in the IQHE. A simple interpretation of the Z2 
invariant is in terms of doublets of edge modes: in the 
QSH phase the edge states consist of an odd number of 
Kramers doublets, whereas the conventional band insula- 
tor is characterized by an even number (including zero) 
of pairs of edge states. The nrlH ^'mbov of Kramm 
doublets is p|jW|St against disorderliZHl3ll3'Ej , ElE3Ej and 
interactions .EZlEa In particular, when subject to time- 
reversal invariant random impurity potentials, there is 
always one perfectly conducting channel, as long as the 
bulk topological properties are not altered by disorder. 

There is a natural generalization oi|-ti|* 
cal insulator to three dimensions (3D).ffiLIJli 

the 2D version, there are now four independent Z2 topo- 
logical invariants which describe the number of Kramers 
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TABLE I: Ten symmetry classes of single particle Hamiltonians classified in terms of the presence or abpapee of time-reversal 
symmetry (TRS) and particle-hole symmetry (PHS), as well as sublattice (or "chiral") symmetry (SLS)£j£j In the table, the 
absence of symmetries is denoted by "0". The presence of these symmetries is denoted either by or "—1", depending 

on whether the (antiunitary) operator implementing the symmetry at the level of the single-particle Hamiltonian squares to 
"+l"or "—1" (see text). [The index ±1 equals r\ c in Eq. (|lb|); here e c = +1, —1 for TRS and PHS, respectively.] For the first 
six entries of the TABLE (which can be realized in non-superconducting systems) TRS = +1 when the SU(2) spin is integer 
[called TRS (even) in the text] and TRS = — 1 when it is a half-integer [called TRS (odd) in the text]. For the last four 
entries, the superconductor "Bogoliubov-de Gennes" (BdG) symmetry classes D, C, Dili, and CI, the Hamiltonian preserves 
SU(2) spin-1/2 rotation symmetry when PHS=— 1 [called PHS (singlet) in the text], while it does not preserve SU(2) when 
PHS=+1 [called PHS (triplet) in the text]. The last three columns list all topologically non-trivial quantum ground states as 
a function of symmetry class and spatial dimension. The symbols Z and Z 2 indicate whether the space of quantum ground 
states is partitioned into topological sectors labeled by an integer or a Z 2 quantity, respectively. 



degenerate band crossings (Dirac points) in the spectrum 
on the surface of the 3D bulk, thereby distinguishing 
the conventional insulator, the topologically trivial phase 
from the topologically non-trivial phase. Although the 
effects of disorder and interactions on the Z 2 topological 
insulator have been less well studied in 3D than in the 2D 
case, there are known to exist gapless surface modes in 
the topologically non-trivial 3D phase which are robust 
against arbitrary strong disorder as long as the latter 
does not alter the bulk topological, pr . pnpr . tiGs , j in analogy 
to the QSH effect (QSHE) in 2DJl3€llESEl]H3£3 These de- 
localized surface states, whose Fermi surface encloses an 
odd number of Dir ac . pum ts, form a two-dimensional "Z2 
topological metal" .E3c3£a 

Recently, a series of experiments have been per- 
formed on certain candidate materials for Z2 topologi- 
cal insulators. For example, the QSH effect has been 
obser ved . in ., HgT c/(Hg,Cd)Te semiconductor quantum 
wells.BBEE3ia Moreover, a 3D Z 2 topological phase 
has been predictetLiap-strained HgTe and for Bismuth- 
Antimony alloys .E3BeJ Indeed, photoemission experi- 
ments on the latter system have revealed an odd number 
of Dirac points inside the Fermi surface on the (111)- 
surface, thereby providing (indirect) evidence for the ex- 
istence of a Jimii-trivial topological phase in three spatial 
dimensions 

In this paper we provide an exhaustive classification of 
topological insulators and superconductors. Our classifi- 
cation is for non-interacting systems of fermions. How- 
ever, since there is a gap, our results also apply to in- 
teracting systems as long as the strength of the interac- 
tions is sufficiently small as compared to the gap. As 
the majority of previous works studied two-dimensional 



topological phases, we shall be mostly concerned with 
the classification of 3D systems, and only briefly com- 
ment on one- and two-dimensi onal t opological insulators 
in the discussion section (Sec. VIII ). In the same spirit 
as in the treatments of Z2 topological insulators, we im- 
pose several discrete symmetries on a family of quantum 
ground states. We then ask if different quantum states 
can be transmuted into each other, without crossing a 
quantum phase transition, by a continuous deformation 
respecting the discrete symmetries. 

If we are to include spatially inhomogencous deforma- 
tions of quantum states, such as those arising, e.g., from 
the presence of random impurity potentials, the natural 
discrete symmetries we should think of would bp-,thosc 
considered in the context of disordered systems £3 It is 
at this stage that we realize that the existence of the 
classification of random Hamiltonians, familiar from the 
theory of random matrices, will become very useful for 
this purpose. 

Specifically [following Zirnbauer, and Altland and Zirn- 
bauer (AZ)£3'EZl all possible symmetry classes of random 
matrices, which can be interpreted as a Hamiltonian of 
some non-interacting fermionic system, can be system- 
atically enumerated: there are ten symmetry classes in 
total. (For a summary, see Table |[) The basic idea as 
to why there are precisely ten is easy to understand. 
Roughly, the only generic symmetries relevant for any 
system are time- reversal symmetry (TRS), and charge 
conjugation or particle-hole symmetry (PHS). Both can 
be represented by antiunitary operators on the Hilbcrt 
space on which the singlc-pa*ticle Hamiltonian (a ma- 
trix) acts, and can be writtenEa on this space in the form 
KU, with K = complex conjugation, and U = unitary. 



Any of these two symmetries can either be absent, which 
we denote by 0, or be present and square to the identity 
operator, or to minus the identity operator, which we de- 
note by +1 or — 1, respectively. This gives nine possible 
choices for the pair of symmetries TRS and PHS (Table 
||). However, since we consider TRS and PHS, we may 
also consider, in addition, their product SLS := TRS x 
PHS, often referred to as "sublatticc" (or "chiral") sym- 
metry. Now, for eight of the nine assignments of a pair 
of 0,±1 to the pair of symmetries TRS and PHS, the 
presence or absence of the product SLS of these symme- 
tries is uniquely determined (Table |). But the assign- 
ment (TRS, PHS) = (0,0) allows for SLS to either be 
present (SLS = 1) or absent (SLS = 0). Therefore one 
obtains ten symmetry classes (Table ||), an exhaustive 
list. 

The so-obtained ten (AZ) symmetry classes of random 
matrices are conventionally named after the mathemati- 
cal classification of symmetric spaces, and called A, AI, 
All, AIII, BDI, CII, D, C, Dili, and CI (Table @). This. 
AZ classification includes the three previously known,c2l 
so-called "Wigner-Dyson symmetry classes" (or "stan- 
dard symmetry classes") relevant for the physics of An- 
derson localization of electrons in disordered solids [corre- 
sponding to orthogonal (AI), unitary (A), and symplectic 
(All) random matrix ensembles] . Three so-called "chiral 
classes" can be obtained, from the Wigner-Dyson classes, 
by imposing an additional SLS; these arc conventionally 
called "chiral orthogonal" (BDI), "chiral unitary" (AIII), 
and "chiral symplectic" (CII) symmetry classes. A well- 
known prototypical example of a system in a chiral sym- 
metry class is a disordered tight-binding model on a bi- 
partite lattice, such as tinp random hopping model, and 
the random flux model.cil Finally, there are four addi- 
tional symmetry classes (D, C, Dili, and CI) describ- 
ing the (Anderson-like) localization physics of the non- 
interacting Bogoliubov-de Gennes (BdG) quasiparticles 
existing deep inside the superconducting state of disor- 
dered superconductors, as described within a mean-field 
treatment of pairing. (Four symmetry classes arise since 
SU(2) spin-rotational invariance, or TRS may be present 
or absent.) 

In terms of this terminology, the above mentioned Z2 
topological insulator is a topologically non-trivial insu- 
lator within the symplectic (or: "spin-orbit") symmetry 
class (class All). In this paper, we pursue this direction 
further and provide a classification of all possible topo- 
logical insulators in 3D. Specifically, we take the clas- 
sification by AZ and ask if different Hamiltonians can 
be continuously deformed into each other within a given 
symmetry class. 



A. Summary of results 

One of the key results of the present paper is our find- 
ing that for five out of the above-mentioned ten (AZ) 
symmetry classes for random matrices, there exist topo- 
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logically non-trivial insulatorsB in three spatial dimen- 
sions. These classes are: 

• symplectic symmetry class (class All), 

• chiral unitary symmetry class (class AIII), 

• chiral symplectic symmetry class (class CII), 

• BdG symmetry class of superconductors with time- 
reversal (TR) but no SU(2) spin rotation symmetry 
(class Dili), 

and 

• BdG symmetry class of superconductors with both 
TR and SU(2) spin rotation symmetries (class CI). 

All these—symmetry classes possess a TRS of some 
form.aoo Our result is summarized in the last col- 
umn of Table [j], where the symbols Z and Z2 indicate, 
whether the space of quantum ground states is parti- 
tioned into topological sectors labeled by an integer or a 
Z2 quantity, respectively. 

We are going to derive these findings by using two 
complementary strategics. First, wc introduce a suit- 
able topological invariant which takes on integer values, 
and can be used to label topologically distinct quan- 
tum ground states (see Sec. |n|). Second, we study two- 
dimensional boundaries terminating 3D topological insu- 
lators and use the appearance of gapless surface modes 
as a diagnostic for the topological nature of the 3D 
bulk properties (see Sec. |TV| ). The latter is accomplished 
by considering Dirac Hamiltonian representatives of 3D 
topological insulators. Before we turn to a more detailed 
and technical discussion, we outline below these strate- 
gies in general terms. 



1. Bulk topological invariant 

To characterize topological properties of bulk wave 
functions in classes AIII, Dili, and CI, which are classes 
that can be realized as time reversal (TR) invariant su- 
perconductors, we introduce an integer-valued topolog- 
ical invariant (winding number), to be denoted by v. 
This winding number can be defined for those symme- 
try classes in which the Hamiltonian can be brought into 
block off-diagonal form, a well known property of ran- 
dom matrices in all the so-called chiral symmetry classes 
AIII, BDI, and CII. This turns out to be also a prop- 
erty of symmetry classes Dili and CI (and arises from 
the presence of a ("sublatticc", or "chiral") symmetry 
which is a combination of PHS and TRS). One of the 
simplest examples of such a quantum state with a non- 
trivial winding numberJsJjji fact the well known Balian- 
Werthamer (BW) stateBEO of the B phase of liquid 3 He. 
In the above language, the BdG fermionic quasiparticles 
in 3 He B are in a 3D topological insulating phase in class 
Dili, with winding number v = 1. 
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However, when additional discrete symmetries arc 
present in a given symmetry class, these can (and will) 
restrict the possible values of the winding number v to a 
subset of the integers. Indeed, while an arbitrary integral 
winding number can be realized in classes Dili and AIII, 
only an even winding number turns out to be allowed for 
class CI. For classes BDI and CII, on the other hand, 
we do not find any example with a non-trivial winding 
number. 



2. Surface Dirac/Majorana fermion modes 

As anticipated from the examples of topological insula- 
tors in 2D such as the IQHE and the Z2 topological insu- 
lators, the non-trivial topological properties of the quan- 
tum state in the 3D bulk manifest themselves through 
the appearance of gapless modes at a 2D surface termi- 
nating the 3D bulk; these turn ■QLyt-jp.l^-ttaaless Dirac 
(or Major ana) fermion modes .OEJCj£JSl3E3 The con- 
verse is also true: the physics at the boundary faithfully 
reflects (indeed " holographically" ) the non-trivial topo- 
logical features of the bulk quantum state, in a way that 
is reminiscpa t of the situation familiar from the quantum 
Hall statesc3'E3E3EZl [and also from recent work on grav- 
itational Chern-Simons theory in (2 + 1) dimensions .L3] 
We are thus led to consider the nature and the properties 
of Dirac fermions appearing at two-dimensional surfaces 
terminating the 3D bulk as a tool to characterize and to 
learn about the topological properties of the wave func- 
tions of the three-dimensional bulk of interest. 

Indeed, a complete classification scheme of properties 
of Dirac fermions in two spatial dimensions has recently, 
appeared in the work of Bernard and LeClair (BL),l£| 
and we will use their results extensivcly.Ea Interestingly, 
and of key importance to our goal, is the fact that the 
BL-classification scheme consists of 13 symmetry classes, 
not just the 10 (AZ) classes mentioned above. This is 
due to the "Dirac structure" of the Hamiltonian: in ad- 
dition to the ordinary ten symmetry classes, each one of 
the three classes AIII, Dili and CI (i.e., classes which 
can be realized in TR invariant superconductors), subdi- 
vides in fact into two symmetry classes (not just one, as 
in the AZ scheme which applies to random Hamiltonian 
not of "Dirac form"). The gapless nature of 2D Dirac 
fermions in the extra three symmetry classes turns out 
to be entirely robust against any perturbation (respect- 
ing the symmetries of a given symmetry class) , including 
those breaking translational invariance (i.e., disorder po- 
tentials). Remarkably, it is precisely these extra three 
classes which are realized at a surface of bulk topological 
insulators in classes AIII, Dili, and CI in three spatial 
dimensions. More specifically, for symmetry classes AIII 
and Dili there exist 3D topological insulators possess- 
ing any number of gapless 2D surface Dirac and Majo- 
rana fermion modes, respectively, which are stable to ar- 
bitrary perturbations (respecting the symmetries). For 
class CI, on the other hand, only an even number of 



such gapless 2D Dirac fermion modes which are robust 
to perturbations can be realized at a sur face .. T his should 
be contrasted with the even-odd cffcct£l't30ll3 governing 
the robustness of the gaplessness of the spectrum of the 
Dirac fermion modes at a surface of a three-dimensional 
Z2 topological insulator, a feature which is protected by 
the Z2 invariant. 

For class CII the situation is similar to the case of 
Z2 topological insulators in class All (the symplectic 
symmetry class). A three-dimensional insulator in class 
CII can support surface Dirac fermions that are stable 
against any symmetry-preserving perturbations. We will 
explicitly demonstrate this below for the case when the 
number of surface Dirac fermions is two. On the other 
hand, when the number of such flavors is twice an even in- 
teger, the surface Dirac fermions are not protected from 
acquiring a mass. This stability of the gapless surface 
Dirac fermions has nothing to do with the winding num- 
ber mentioned above. 



3. Examples and many-body wave functions 

The discussion so far has been solely in terms of single- 
particle physics. This is not to say, however, that these 
results are completely irrelevant to interacting many- 
body systems: in particular, when viewed as mean-field 
ground states, the families of states considered above can 
naturally arise as a consequence of strong correlations. 
For example, the most interesting topological features 
of the Moore- Read Pfaffian state of the fractional quan- 
tum Hall (FQH) effect for the half-filled Landau level 
can be understood in terms of the ground state of a 
certain (p x + ip y ) BCS superconductor, within a mean 
field treatment of pairing £j Moreover, once the dynam- 
ical fluctuations of the pairing potential and of the elec- 
tromagnetic U(l) gauge field are included, the suoeceotu 
ducting ground state is topologically ordered. E3Ej|£jLJ 
In a similar fashion, the 3D topological phases realized 
in superconducting classes Dili, AIII and CI are, when 
viewed as many-body wave functions, 3D analogues p£ 
paired FQH states, such as the Moore-Read PfaffiaH,L3 
the Halperin 33lJ^3 and the Haldane-Rezayi states£3 In 
particular, the BW state of the B-phase of liquid 3 He 
can be thought of as a 3D analogue of the Moore-Read 
Pfaffian state. Below, we will discuss in more detail the 
properties of real-space wave functions of topologically 
non-trivial and of topologically trivial phases of these 3D 
topological insulators, which arc analogues of the familiar 
"weak-pairing" and the "strong-pairing" states, respec- 
tively, of the Moore-Read Pfaffian state in two spatial 
dimensions. We will also discuss topological degenera- 
cies arising in such 3D phases. 
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B. Outline 



This paper is structured as follows. Since the symme- 
try classificationEZrO of Hamiltonians which emerged in 
the context of random matrix theory is crucial for our dis- 
cussion, we give a rather pedagogical description of it in 
Sec. |n]. The topological winding number v is introduced 
in Sec. |l| to characterize the bulk ground state wave 
functions. We have a close look at the surface Dirac. 

following Bernard and LeClair.E2l 



[V 



fermions in Sec. 
Explicit examples of 3D Dirac insulators with a mini- 
mal number of surface Dirac fermions arc constructed in 
Sec. [v|. To clarify the connection between the topological 
properties of the quantum state in the 3D bulk and the 
appearance of gapless Dirac fermions at 2D surfaces, we 
consider a topological field theory description in terms 
of (doubled) Chcrn-Simons theory in Sec. VI. Finally, we 
discuss the 3 D to pological insulators as many-body sys- 
tems in Sec. |VIl| and study their many-body wave func- 
tions and topological degeneracies. We conclude in Sec. 
VIII with a brief discussion of the close connection be- 



tween the non-trivial topological characteristics of the 3D 
bulk and the Anderson (de)localization physics at the 2D 
boundary. 

For those readers who arc interested in details and pre- 
fer a systematic presentation, we recommend to read all 
sections sequentially. On the other hand, those readers 
who prefer to understand the concepts rather through 
explicit examples, may skip Sees. W and fVl and should 
proceed to read only Sees. |ll4 jvj. 



VII 



and VIII. 



II. SYMMETRY CLASSIFICATION OF 
NON-INTERACTING HAMILTONIANS 

We start by recalling-thc basic ideas underlying the 
classification by AZEDO of non-interacting fermionic 
Hamiltonians ( "random matrices" ) in terms of ten sym- 
metry classes. These are classified in terms of the pres- 
ence or absence of certain discrete symmetries (see Table 
H for a summary). The corresponding symmetry opera- 
tions are classified into two types, 

P: H = -PHP- 1 , PP f = l, P 2 = 1, (la) 
C : H = e c CH T C- 1 , C& = 1, C T = i h C, (lb) 

where TL is a matrix or an operator representing a single- 
particle Hamiltonian, e c = ±1 and i] c = ±1. 

The symmetry operation corresponding to C (C-type 
symmetry) represents a TRS operation when e c = 1 and 
a PHS symmetry operation when e c = — 1. Further- 
more, we distinguish two cases, r\ c = ±1: (e c ,?/ c ) = (1, 1) 
represents a TRS for spinless (or integer spin) particles, 
whereas (e c ,?7c) = (1> — 1) represents a TRS for spinful, 
half-integer spin particles. Similarly, (e c ,7y c ) = (—1,1) 
represents a PHS for a triplet pairing BdG Hamiltonian 
whereas (e c ,r] c ) = (—1,-1) represents a PHS for a sin- 
glet pairing BdG Hamiltonian. Note that although the 



form of a C-type symmetry operation can be changed 
by a unitary transformation, the value of (e c , rj c ) remains 
unchanged. The presence of TRS for half-integer spin im- 
plies Kramers degeneracy whereas the presence of PHS 
implies that the energy spectrum is symmetric about zero 
energy. 

Similar to PHS (C-type symmetry with e c = —1), a 
P-type symmetry implies a symmetry of the energy spec- 
trum. In condensed matter systems, it is often realized as 
a sublattice symmetry on a bipartite lattice (i.e., the sym- 
metry operation that changes the sign of wave functions 
on all sites of one of the two sublattices of the bipartite 
lattice) and is sometimes called chiral symmetry. When 
we have two P-type symmetries, P and P', say, we can 
construct a conserved quantity by combining them, i.e. 
[H, PP'} = where PP'(PP')^ = 1. Note that the latter 
property implies that we can block-diagonalize Ti and ap- 
ply our classification scheme to each of the blocks. Con- 
sequently, it is enough to consider ensembles of Hamilto- 
nians possessing only a single, or no P-type symmetry. 

Note that whenever a Hamiltonian possesses both P- 
and C-type symmetries, it automatically has another, 
different C-type symmetry C defined by 



H = e'C'H T C' 



C" = PC, e' = 



(2) 



Thus, symmetry classes of Hamiltonians possessing both 
P- and C-type symmetries automatically possess in fact 
all three, chiral, particle-hole and time-reversal symme- 
tries. 

The complete classification in terms of the presence or 
absence of chiral, particle-hole, and time-reversal symme- 
tries is summarized in Table The ten symmetry classes 
of AZ can be grouped into three categories: the three 
standard (or "Wigner-Dyson" ) classes {A, AI, All}, the 
three chiral classes {AIII, BDI, CII}, and the four BdG 
(superconductor) classes {D, C, DHL CI}. We wish 
to point out, however, that classes CI and Dili can be 
thought of as a close cousin of the chiral classes, since in 
each of these two classes one can find a unitary matrix, 
by combining the TRS and PHS, which anticommutes 
with all members of the class. Conversely, class AIII can. 
also be thought of as a BdG (superconductor) class .cl 
(This will be discussed further below). 

Below, we will give a more detailed and physical de- 
scription for each class. In this section, we use four sets 
of standard Pauli matrices s M , c M , i M , and r M (where 
fi = 0,x,y,z, and sq, cq, to, and tq are 2x2 unit ma- 
trices.) Unless otherwise specified, the Pauli matrices 
act on spin indices (T/I)j whereas c M act on the two (vi- 
and B-) sublattice indices of a bipartite lattice; the Pauli 
matrices are used to represent the particle-hole space 
appearing in the BdG Hamiltonian for quasiparticles in 
a superconductor, whereas are used for superconduc- 
tors for which the 2-component S z of spin is conserved. 
We also use two additional sets of Pauli matrices, er^, 
Ty, (where /i = 0,x, y,z, and <jq and To are 2x2 unit 
matrices). 
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A. Standard (Wigner-Dyson) classes 

Let us first review thu. familiar, standard Wigner- 
Dyson symmetry classes O An ensemble of Hamiltoni- 
ans without any constraint other than being Hermitian 
is called the unitary symmetry class (class A). Imposing 
TRS for half-integer spin on the unitary symmetry class, 

is y H T (—is y ) = 7i, (3) 

one obtains the symplectic symmetry class (class All); 
note that is y is antisymmetric, (is y ) T = —is y . Impos- 
ing, in addition, SU(2) spin rotation symmetry on the 
symplectic symmetry class, one obtains the orthogonal 
symmetry class (class AI), 

H T = H. (4) 

The symmetry operation in Eq. (|]) is the TRS operation 
for integer-spin or spinless particles. We distinguish the 
two C-type symmetries in Eqs. (||) and (H) by referring 
to them as "TRS (odd)" and "TRS (even)" , respectively. 



B. Chiral classes 

Symmetry classes of Hamiltonians possessing a P-type 
symmetry implemented by a unitary transformation 

c z Hc z = -H, (5) 

[letting c z := P in Eq. ([!])] are conventionally called chi- 
ral classes. As already mentioned above, the correspond- 
ing unitary transformation is typically implemented as a 
sublattice symmetry on a bipartite lattice. Equations (|l|) 
and (||) imply that all the energy eigenvalues appear in 
pairs (with a possible exception at zero energy). From 
an eigenstate tjj with energy E, one can obtain another 
state with the opposite energy — E by a unitary transfor- 
mation, C z l\). 

In complete analogy with the standard ("Wigner- 
Dyson") classes discussed above, the ensemble of chiral 
Hamiltonians without any further conditions is called the 
chiral unitary class (class AIII). Imposing TRS for half- 
integer spin (^), we obtain the chiral symplectic class 
(class CII). Imposing both the TRS and SU(2) symme- 
tries [i.e., imposing the TRS (Q)], we obtain the chiral 
orthogonal symmetry class (class BDI). 

A well-known physical realization of the chiral unitary 
symmetry class (class AIII) is a disordered tight-binding 
model on a bipartite lattice with broken TRS, such as the 
random flux problem. EJ However, a class AIII Hamilto- 
nian cancJ also be interpreted as an ensemble of BdG 
Hamiltonians that have a TRS and are invariant under a 
U(l) subgroup of SU(2) spin rotation symmetry (rotation 
around z-componcnt of spin, say), as will be discussed 
further below. 



C. BdG classes 

Following Altland and ZirnbauerJE] we now consider 
the following general form of a Bogoliubov-de Gcnncs 
Hamiltonian for the dynamics of quasiparticles deep in- 
side the superconducting state of a superconductor 

H = \ ( ct, c ) Ha (^ t ) . W 4 = ( _%r ) , (6) 

where H4 is a AN x AN matrix for a system with N 
orbitals (lattice sites), and c = (cf , cjj. [c and can be 
either column or row vector depending on the context.] 
Because of S = (hcrmiticity) and A = — A T (Fermi 
statistics), the BdG Hamiltonian (||) satisfies 

(a) : H^ = -t x Hjt x , [PHS (triplet)]. (7) 

This is a C-type symmetry with (e c , r) c ) = (—1, +1), and 
will be called the PHS (triplet). 

In terms of the presence or absence of TRS (odd) rep- 
resented by 

(b) : H 4 = is y Hj(-is y ), [TRS (odd)], (8) 
and of SU(2) spin rotation symmetry represented by 

(c) : [7U, Ja] = 0, Ja := ( °° _° g T ) , 

a = x,y,z, [SU(2) symmetry] , (9) 

BdG Hamiltonians (Q) are classified into the four sub 
classes listed in Table |: Classes C and CI are primar- 
ily relevant to singlet SC whereas classes D and Dili are 
primarily relevant to triplet SC, although one can also 
consider admixture of singlet and triplet order parame- 
ters in the absence of the parity symmetry, as known in, 
e.g., CcPt 3 Si.E3 

1. BdG classes without spin rotation symmetry 

Class D First, we consider the symmetry class with 
neither TRS nor SU(2) invariance. In this case, a set of 
BdG Hamiltonians satisfying (a) is nothing but the Lie 
algebra so(4m). Any element H.4 £ so(4m) can be diago- 
nalizcd by an SO (Am) matrix g as gH^g -1 — diag(e, —s), 
with e = diag(ei,£2, • • • ) ; i-e., the spectrum is particle- 
hole symmetric. 

An example of class D BdG Hamiltonian in 2D is a 
2D spinless chiral p-wave (p ± ip-wave) superconductor, 
which can be written in momentum space as 

"-5E(4. =-*)"«(£)■ 

h(k)=A(k x t x + k y t y )+e k t z , (10) 

where k = (k Xl k y ) is the 2D momentum, A 6 R is the 
amplitude of the order parameter, and denotes the 
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AZ class 


SU(2) 


TRS 


Constraints on Hamiltonians 


Examples in 2D 


D 


X 


X 


txH. 1 t x = —H 


Spinless chiral (p ± jp)-wave 


Dili 


X 


o 


t x n r t x = -n, isyH 1 ' {-isy) = n 


Superposition of (p + ip)- and (p — ip)-wave 


A 


A 


X 


no constraint 


Spinful chiral (p ± ip)-wave 


AIII 


A 


o 


TyHr-y = — H 


Spinful p x - or p^-wave 


C 


O 


X 


r y n r r y = -n 


(d ± id)-wave 


CI 


O 


o 


r y H r r y = -H, H* =H 


d x 2_ y 2- or d xy -wa,ve 



TABLE II: Symmetry classification of Bogoliubov-de Gennes systems, in terms of the presence ("O") or absence ("x") of 
SU(2) spin rotation symmetry and time-reversal symmetry (TRS). In classes A and AIII. Hamiltonians are invariant under 
rotations about the z- (or any fixed) axis in spin space, but not under full SU(2) rotations, as denoted by "A" in the table. 
The sets of standard Pauli matrices s x ,y,z, t x ,y,z, and r IjS , z act on the spin, the particle-hole, and the grading defined in Eq. 
(Eq), respectively. 



energy dispersion of a single particle. The Hamiltonian 
has PHS (triplet) [Eq. (@)], h(k) = -t x h T {-k)t x . 

Class Dili Consider class Dili, which satisfies con- 
ditions (a) and (b). A set of matrices which simultane- 
ously satisfy (a) and (b) does not form a subalgebra of 
so(4m), but consists of all those elements of the Lie alge- 
bra so(Am) which are not elements of the sub Lie algebra 
u(2to).EI 

Combining (a) and (6), one can see that a member of 
class Dili anticommutes with the unitary matrix t x ®s yi 



H± = —t x <8> SyHitz 



(11) 



It is interesting to note, that Hamiltonian (fill) is 
a direct product of Hamiltonian (|l0|), h(k x ,k y )~^ and 
h(—k x ,ky). This follows from a simple reordering of 
the basis elements in Eq. (|l3|), such that (cj c ,c_ fe ) — > 
( c fct! c -fcT> c ki> c -fcj.)- The superconductor described by 



the order parameter Eq. (13b) can be thought of as 
a two-dimensional analogue of the BW state realized 
in the B phase of 3 He. The BW state, which is also 
a member of class Dili and described by the d-vector 
dk = A ( k x , k y , k z ) , will be discussed in Sec. [V] as an 
example of 3D topological insulators. 



In this sense, class Dili Hamiltonians have a chiral struc- 
ture. It is sometimes convenient to take a basis in which 
the chiral transformation, which is t x ® s y in the present 
basis, is diagonal. In one of such bases, a class Dili 
Hamiltonian takes on the form 



Ha 



D 
£>t 



D 



-D 



(12) 



An example of a 2D BdG Hamiltonian in symmetry 
class Dili is a p x -wave (or p y -wave) superconductor with 
e?-vcctor not pointing the z direction. An equal superpo- 
sition of two chiral p-wave SCs.wJt.h opposite chiralities 
(p x + ip y - and p x — zp y -wave)oOOO also falls into this 
class. The latter can be explicitly written in momentum 
space as 



2. BdG classes with S z conservation 

Let us consider BdG Hamiltonians which are invariant 
under rotations about the z- (or any fixed) axis in spin 
space, yielding to the condition \Ha, Jz] = 0, which im- 
plies that the Hamiltonian can be brought into the form 



Ha 



( 


a 








b 







a' 


-b T 










-b* 







\ 










-a lT 



J 



(14) 



Due to the sparse structure of H.4 , we can rearrange the 
elements of this 47V x 4N matrix into the form of a 2A^ x 
2A^ matrix, 



,(13a) 



with the row vector (c\,c_ k ) = (cj^, cj^, c_ k p c_ kl ), and 



the matrix elements 



^ / ^ iky 

■ ~\~ ^kf. 



A fe = (dk ■ s)(is y ). 



(13b) 



In terms of the d- vector, dk — A ( — k x , k v , ), the su- 
perconducting order parameter reads 



(13c) 



H 



a 
6t 



-,'T 



1 r / 

+ -tr W 



(15) 

Note that the Hamiltonian H is traceless. 

To summarize, a BdG Hamiltonian which is invariant 
under rotations about the z-component in spin space can 
be brought (up to a term which is proportional to the 
identity matrix) into the form 



H 



( C T> c i) 



Ho 



Ho 



6 



^ l , (16) 



where £j = Without further constraints, this Hamil- 
tonian is a member of class A (unitary symmetry class). 
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A physical realization of this class is a 2D spinfull chiral 
p-wave (p ± ip-wave) superconductor [c.f. Eq. (|l3|)l with 
d-vector parallel to the z-direction. 

d k = zA(k x + ik. y ) = A ( 0, 0, k x +ik y ) . (17) 

More compactly, the spinful chiral (p ± ip)-wavc super- 
conductor can be expressed as 



H = 



where 



E 



-hi 



4 ~tl- 



'fc1 
t 

-U 



S k = A(k x + iky), 



,(18) 



(19) 



and £f/l,k = ^T/i-fc' S k = S -k ( b y definition). [Inci- 
dentally, the single-particle Hamiltonian defined by Eqs. 
( |l8| ) and (|l9| ) happens to belong to class D, if ^ k = 
The Anderson- Brinkman-Morel (ABM) state, which is 
3D superfluid given by the same order parameter as 2D 
spinful chiral p-wave superconductor [Eq. (17)], is also a 
member of class A. 

Class C (full SU(2) symmetry) If we further im- 
pose the full SU(2) rotation symmetry, £ CT and d are con- 
strained by 

£i = £t=:£, 5 = S T . (20) 

These two conditions can be summarized as 

r v n^r v = -H 2 , [PHS (singlet)]. (21) 

This is a C-type symmetry with (e c , r/ c ) = (—1,-1), and 
will be called the PHS (singlet). 

An example of a BdG Hamiltonian in class C ki-two 
dimensions is the (d + id)-w&ve superconductor.LJ Its 
Hamiltonian is given by Eq. (|l8|) together with the ma- 
trix elements = ^ k = £fc and 



^^xy^x^y i 



(22) 



where A x 2_ y 2 and A xy are real amplitudes for the d x i_ y i 
and id xy superconducting order parameters, respectively. 

Class CI (full SU(2) symmetry + TRS) Imposing 
both full SU(2) rotation and TR symmetries, leads, in 
addition to the constraints (|20|), to 



C = £, 6* = 5. 
These conditions can be summarized by 



r y H 2 r y = -n 2 



Ti 9 = 



(23) 



(24) 



Combining these two conditions we can obtain a P- 
type (i.e., chiral) symmetry, r y Ti.2r y = —Jii. It is also 
convenient to rewrite the Hamiltonian by rotating the r M 
matrices by (r x ,r y ,r z ) => {r x ,—r z ,r y ). In this basis, tJbe 
class CI Hamiltonian takes on block off-diagonal formES 



n 2 = 



D 
Dt 



D = S - i£ = D 



(25) 



An example of a BdG Hamiltonian in class CI in two 
dimensions is a 2D d-wave (d 2 ,2_ J/ 2-wave) superconductor, 
which is described by the order parameter Eq. (123), with 



A 



xy 



0, but with A x 2_ y 2 ^= 0. 



Class AIII (S z conservation + TRS) If, in addi- 
tion to conservation of the ^-component of spin we fur- 
ther impose TRS (b), we obtain the conditions 



which can be summarized as 



r v H 2 r v 



5 = 8\ 



-Hi 



(26) 



(27) 



By interpreting the r-grading as c-grading, this is noth- 
ing but the SLS. Thus, a BdG Hamiltonian possessing 
TRS (b) and conserving one component of SU(2) spin, 
can be thought ,-ef as a member of the chiral class with- 
out TRS (AIII)J±a An example of a BdG Hamiltonian of a 
two-dimensional superconductor in class AIII is a p-wave 
(Px,2/-wave) superconductor with the d-vector parallel to 
the ^-di recti on, i.e., Hamiltonian (|l3] ) with d k = Azk. 



in Eq. (13c), or, alternatively, Hamiltonian ( |18[ ) with 

Sk = Ak x ,y 

Before closing this section, we emphasize that the sym- 
metry classification we have described also applies to in- 
teracting fermion systems, as discrete symmetries can be 
imposed on the second quantized fermion creation and 
annihilation operators. 



III. CHARACTERIZATION IN THE BULK 

A useful quantity to discuss the bulk characteristics 
of topological insulators is the spectral projection oper- 
ator. We will first discuss the spectral projector for a 
general Bloch Hamiltonian with a bulk band gap, and 
then specialize to those symmetry classes that satisfy a 
certain discrete unitary symmetry (called P-typc in Sec. 
||), in which case the projection operator can be brought 
into block off-diagonal form. This block off-diagonal rep- 
resentation of the projector allows for the definition of 
a winding number that distinguishes between different 
topological phases. 



A. Projection operator 

In the presence of translation invariance, ground states 
of non-interacting fermion systems can be constructed 
as a filled Fermi sea in the d-dimensional Brillouin zone 
(BZ), in Fourier space. The band structure can then be 
viewed as a map from the BZ to the space of Bloch Hamil- 
tonians. Similarly, the spectral projection operatoiO can 
be thought of as a map from the reciprocal unit cell to 
a certain Lie group or cosct manifold, which wc will call 
space of projectors or target space. In order to define the 
spectral projector, let us consider an eigenvalue problem 
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at momentum k, 



H(k)\Ua(k))=Ea(k)\Ua{k)), 



(28) 



where TL{k) is a Hamiltonian at k in the BZ, \ua{k)} is 
the Bloch state in the a-th band with energy Ea(k). We 
assume the existence of a bulk band gap centered around 
some energy Eq and define the quantum ground state by 
filling all states with Ea(k) < Eq . Without loss of gener- 
ality, we can set E = 0, as a suitable constant chemical 
potential in the definition of single-particle energy lev- 
els can always be included. The projector onto the filled 
Bloch states at fixed k is then defined as 



filled 



p(k) = J2 M*))M*)|. 



It is convenient to introduce 



Q{k) = 2P(k) - 1. 



(29) 



(30) 



It is readily checked that the so-defined "Q-matrix" is 
characterized by the conditions 



= Q, Q 2 = 1, trQ = m — n, 



(31) 



where we consider the situation where, at each fc, we 
have m filled and n empty Bloch states. Depending on 
the symmetry class, additional conditions may be im- 
posed on Q, which we will consider later. Without any 
such further conditions, the projector takes values in the 
so-called Grassmannian G m . m +„(C): the set of eigenvec- 
tors can be thought of as a unitary matrix, a member 
of U(m + n). Once we consider a projection onto the 
occupied states, we have a gauge symmetry U(m) for 
the occupied states, and a similar gauge symmetry U(n) 
for the empty states. Thus, each projector is described 
by an element of the coset U(to + n)/U(m) x U(n) ~ 
G m ,m+n{'£) — G nim + n (C). On the other hand, an ele- 
ment of G m ^ m + n (<C) can be written as 

Q = UKU\ A = diag(I m , -I„), U £ TJ(m + n). (32) 

(where m eigenvalues of the diagonal matrix A equal +1, 
and the remaining n eigenvalues equal — 1.) Imposing 
TRS or PHS, which are realized by an anti-unitary op- 
eration (called C-type in Sec. |l|), prohibits certain types 
of maps from the BZ to the space of projectors (see Ta- 
ble |nij). On the other hand, if the discrete symmetry 
is realized by a unitary operation (called P-type in Sec. 
H) the space of projectors is altered from G n , m + n (C) to 
U(m) (see the following subsection). 

We now ask if any clement of the set of projec- 
tors within a given symmetry class can be continu- 
ously deformed into any other, without closing the en- 
ergy gap. Mathematically, this is related to the homo- 
topy group of the (topological) space of projectors. For 
two spatial dimensions the relevant homotopy group is 
^2[G m ,m+n{ < C)] = implying that the projectors are 



classified by an integer (Chern number): IBB projectors 
with different Chern numbers cannot be deformed into 
each other adiabatically; this is the mathematical reason 
why there exists a series of distinct 2D quantum Hall 
insulators labeled by an integer, which is the Hall con- 
ductivity a xy (measured in natural units of e 2 /h). On 
the other hand, in three spatial dimensions and in the 
absence of additional— discrete symmetries, the relevant 
homotopy group isLS'O 



n 3 [G 



m,m+n 



(C)] * {e}, 



(33) 



where {e} represents a group with only one element 
(identity \_ Thus, there is no notion of winding in 3D in 
this casc.El [When m = 1, there is an accidental winding, 
^[Gi, 2 (C)] = Z, because of the Hopf map.BEI] 

This is not to say that there is no topological distinc- 
tion when additional discrete symmetries are imposed on 
the projector. For example, in class All (the symplec- 
tic class), the projector must also satisfy the condition 
(i<j y )Q{k) T (— i<jy) = Q(—k), which arises from the pres- 
ence of TRS. Due to this additional constraint, two dif- 
ferent Q-Held configurations might not be continuously 
dcformablc into one another. Properties of the projector 



for each class are summarized in Table ED . 



B. Block off-diagonal projection operators 

A symmetry realized by an anti- unitary operation (i.e., 
TRS or PHS) relates the projector at wavevector k and 
the one at wavevector —k. Thus, the role of PHS or TRS 
is to prohibit certain Q-ficld configurations in momentum 
space, or to change the topology of the BZ by orbifolding 
k = -k. 

In contrast, imposition of a discrete symmetry which 
is realized by a unitary operation (i.e., SLS or a product 
of PHS and TRS) imposes a condition on the projector 
at each k. Thus, the role of this type of symmetry is to 
change the target manifold of the projector. The main 
focus below is on the projector in those symmetry classes 
which possess (in some basis) a block off-diagonal repre- 
sentation of the Hamiltonian and of the projector, due to 
the presence of a chiral symmetry. This is the case for 
all three chiral classes, AIII, BDI, and CII, and for two 
of the four BdG classes, classes CI and DHL For these 
symmetry classes, the Q-matrix can be brought (upon 
basis change) into block off-diagonal form 



Q 



q 
gt 



(34) 



Since Q 2 = 1, one has qq^ = q^q = 1, and thus q is a 
member of U(m). As before, one should bear in mind 
that the g-matrix can further be subject to several addi- 
tional constraints coming from additional discrete sym- 
metries imposed on Hamiltonians. Properties of the pro- 



jector for each class arc summarized in Table III 
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BUCILL Ul U1U KyLlUlo 111 llHJlllt,llL U.111 BUOL.C 


BL 




fp vim/in i r* rpn ips 

llllUlllLj 1 UllL-Cl 




class 




class 




NLaM target space 


WZW term 


A 




n 


i 
i 


TTfOTVVTTY /V"! v TIY/Vl 


r I UloKcil 


AI 


i O(k) G G,„ mJ _^ fCI 1 OCfcl* — Of-fc) > 


4_i_ 


2 


So^AO/SofAM x Su(N) 


N/A 


All 


\Q(k) G G 2m 2(m+n)(C) 1 (io-^lQffc)*^-^,,) = Q(-fc) j 


3+ 


1 


0(2N)/0(N) x O(iV) 


Z 2 


A TTT 


j <7(/cJ € V(m) } 


± Or z 


i or z 




VV Zj vv 


BDI 




Q i 


2 




N /A 


CII 


i a(k) G U(2m) 1 ficr„Wffc)*f-zcr„) = a(-k) \ 


9 


2 


U(2JV)/0(2AT) 


Z 2 


D 


{Q(fc) G G m , 2m (C)|r a; Q(A;)*r !C = -Q(-fc)} 


3 


1 


0(2Af)/U(A) 


Pruisken 


C 


{ Q(k) e G m , 2m (C) | TyQ(k)*T y = -Q{-k) } 


4_ 


2 


Sp(JV)/U(AT) 


Pruisken 


Dili 


{q(k)eV(2m)\q(ky = -«(-*)} 


5 or 7 


1 or 2 


0(2N) x 0(2AT)/0(2iV) 


WZW 


CI 


{?(fc) GU(m)|g(fc) J = g(-fc)} 


6 or 8 


2 or 4 


Sp(JV) x Sp(AT)/Sp(JV) 


WZW 



TABLE III: The space of projectors in momentum space for each Altland-Ziwibauer (AZ) class. The BL classes represent 
the classification of 2D Dime Hamiltonians obtained by Bernard and LeClairJ£3 and A r ™ m is the smallest possible number of 
flavors of 2D two-component Dirac fermions. The fermionic replica NLaM target spaces, with possible 2D critical behavior [in 
terms of whether it is possible for a given NLcrM to have a topological term, either of Pruisken (IQHE) type or Z 2 type, or a 
Wess-Zumino-Witten (WZW) term] are also listed according to Refs. g and 



The relevant hornpiocy group for projectors that take 
the form Eq. © isBEJ 

!{e}, for d even, 
(35) 
Z, for d odd, 

for m > (d+ l)/2, instead of 7r d [G m , Tn+ „(C)]. 

If a symmetry class of interest does allow topologically 
non-trivial configuration of the projector, a useful tool to 
investigate if a given quantum ground state belongs to a 
non-trivial topological sector (in the space of projectors), 
is a quantized invariant. In the IQHE, it is the quan- 
tized Hall conductivity a xy , which is essentially a winding 
number characteristic of 7r2[G 7Ti . m +„(C)] = ^nfthl ^fpj fe 
topological insulators, it is the Z2 invariant ,EIBlIOI1j'l2I 
which can be constructed from the SU(2) Wilson loops 
that are quantized because of the TRS. Below, we will 
introduce a topological invariant that is applicable for 
symmetry classes with block off-diagonal projector. 

C. Winding number in 3D 

Since 7T3[U(to)] ~ Z (for m > 2), there are, in three 
spatial dimensions, topologically non-trivial configura- 
tions in the space of projectors for those symmetry classes 
for which the Q-matrix can be brought into block off- 
diagonal form. To characterize these distinct classes, we 
can define the winding number 

/d 3 k 

(36) 

where q(k) G U(m), fi,v,p= k x ,k y ,k z , and the integral 
extends over the entire Brillouin zone for lattice systems, 
which is the three-torus T 3 , whereas for continuum mod- 
els the domain of integration in Eq. ([36]) is topologically 
equivalent to three-sphere S* 3 . 



In class AIII, the winding number v[q] can be any inte- 
ger. Due to additional constraints on the g-field in classes 
Dili, CI, BDI and CII, all integers might not be realized. 
One way to determine which integer values can be real- 
ized as winding numbers in each symmetry class, is to 
count the number of flavors of massless Dirac fermions 
allowed by the symmetries at the surface. This counting 
will be done in the next section, yielding the following 
results. 

An arbitrary number of flavors can be realized in sym- 
metry classes AIII and DHL The gaplcss nature of the 
surface Dirac (Majorana) fermions is, irrespective of the 
number of flavors, stable against any perturbations re- 
specting the sjpametries of a given class, including disor- 
der potentials.^ 3 ] 

On the other hand, only an even number of flavors is 
allowed in classes CI, BDI and CII. This suggests that 
the winding number for classes CI, BDI and CII can take 
on only even integer values. Indeed, in Sec. 0, we will 
construct an explicit example of a topological insulator in 
class CI in three spatial dimensions, with winding num- 
ber v = 2 and with two flavors of surface Dirac fermions. 

While an even number of flavors of surface Dirac 
fermions might appear, at first sight, to be possible in 
class BDI, a detailed study of the form of generic pertur- 
bations reveals that the gapless nature of Dirac fermions 
at the surface is not protected in this symmetry class. 
This suggests that the space of all Hamiltonians in class 
BDI has no non-trivial topology. [See also the discussion 
around Eq. @.] 

The stability against perturbations of the gaplcss 2D 
surface Dirac fermions in class CII depends on whether 
the number of flavors is an even or an odd multiple of 
two: if it is even (odd) multiple, the gapless spectrum is 
unstable (stable). We thus expect that the space of all 
Hamiltonians in class CII has a Z2 classification, as in 
class All (the symplectic symmetry class). Furthermore, 
this Z2 classification has nothing to do with the winding 
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number, as we will demonstrate in Sec. ^ where we con- 
struct a 3D class CII insulator with two flavors of surface 
Dirac fcrmions yet with vanishing winding number. 

Before closing this section, several comments are in 
order, (i) Similar discussions are possible in all odd spa- 
tial dimensions, where a winding number can always be 
defined for classes AIII, BDI, CII, CI and DHL In par- 
ticular, since 7Ti [U(m)] ~ Z, topologically non-trivial in- 
sulators characterized by an integer invariant can exist 
in one spatial dimension (ID) when there iSj-a. SLS.L.U 
On the otbp, hand, in the presence of PHS,EjQE!3 or 
TRS (odd), El the non-trivial topological features in the 
bulk of ID topological insulators can be characterized by 
U(l) or SU(2) Wilson loops of the Berry connection, re- 
spectively. Although the quantized values of the Wilson 
loops are not affected by the choice of Bloch wave func- 
tions, as different Bloch wave functions arc related by a 
gauge transformation, the quantized values of the WiL. 
son loops do depend on the choice of the unit cellUj'LlI 
When there is translation invariance, the choice of unit 
cell is arbitrary, whereas if we introduce a boundary, the 
choice of unit cell should be consistent with the loca- 
tion of the boundary. It is in this sense, that the (quan- 
tized) values of Wilson loops in ID reflect the bound- 
ary physics, and are not solely determined from the bulk 
properties, (ii) The relevant homotopy group governing 
the existence of topological insulators in two spatial di- 
mensions is 7T2[U(m)] = {e}. This immediately tells us 
that there are no topological insulators in class AIII in 
2D. However, due to constraints arising from the presence 
of additional discrete symmetries there is still the possi- 
bility of having 2D topological phases in other "crural" 
classes. Indeed, in class Dili, for example, one can con- 
struct a topological insulator ifrop. the mixture of p + ip 
and p — ip pairing states.ESII3O'0|This state is a direct 
analogue of the Kane-Mele modelQ on the honeycomb 
lattice, which-js the mixture of the two Haldane models 
of the IQHE.E3 (iii) Finally, the winding number defined 
above in momentum space (which requires translational 
invariance) can also be defined for disordered systems, 
in a similar fashion in which the Chern number can be 
defined for disordered systems SB 



IV. CHARACTERIZATION AT THE 
BOUNDARY 

A physical consequence of the non-trivial topological 
properties of the quantum state in the bulk is the ap- 
pearance of the gaplcss boundary modes. (Some explicit 
examples will be constructed in the next section.) Con- 
versely, most of the possible bulk phases in (3+1) dimen- 
sions can be inferred by studying their possible (2+1)- 
dimensional boundary physics. In thia-section we con- 
sider, following Bernard and LeClairpJ the symmetry 
classification of 2D Dirac Hamiltonians of the form 



where V± = V± and = A_ and d = d x — id y , 8 = 
d x +idy. Possible dimensionalities of the matrices V± and 
A± depend on the symmetry class as we will see below. 

As before, we impose two types of discrete symmetries, 
P and C. The form of the matrices P and C, acting by 
conjugation, is constrained by the requirement that they 
do not change the kinetic term in Eq. (|37j), resulting in 
the following block diagonal form 



P = 



7 

-7 



C = 



a 
-e c <7 



(38) 



where 7 and a are a matrix satisfying 

77* = 1) 7 2 = 1) cr(jt = 1) ° T = -Vc£ c v- (39) 

All possible forms of 7 and a are listed in Ref. |5J]. 
Due to the Dirac kinetic term, the Bernard-LeClair (BL) 
classification is finer than the ten symmetry classes of AZ: 
there are 13 symmetry classes denoted by 0, 1, 2, 3±, 4±, 
5, 6, 7, 8 and 9±. While all AZ symmetry classes except 
AIII, Dili and CI are in one-to-one correspondence with 
the BL classes, two BL classes correspond to each of the 
AZ symmetry classes AIII, Dili and CI: BL classes 1 and 
2 correspond to AIII, BL classes 5 and 7 to Dili, whereas 
BL classes 6 and 8 correspond to CI. (For a summary, 



see Table [II 



Of direct relevance to our discussion of 3D topologi- 
cal insulators is the minimal number iV™ m of flavors in 
the BL classification. For BL class 3. 



AZ classification), N™ ln is 



H = 



V+ + V- -id + A, 
-id + A_ V+- V- 



(37) 



_l_ (class All in the 
j id ui±c. Since a single flavor of 
2D gapless Dirac fermion cannot be realized on a 2D 
lattice without breaking TRS, the case with Nf in = 1 
(or NJ 1 " 1 = odd) in class 3+ should correspond to a 
state appearing at the two-dimensional surface of a three- 
dimensional Z2 topological insulator in class AIL This 
situaUoH^is indeed realized in the model of Fu-Kanc- 
Mele.O'tHl For BL classes 1, 5 and 6, the minimal number 
TV™ 111 of Dirac fermion flavors is half the minimal number 
of flavors required for the BL classes 2, 7 and 8, respec- 
tively. (Compare e.g. Table [II.) We thus expect that 
2D Dirac fermion modes in the BL classes 1, 5 and 6 are 
realized as boundary states of a non-trivial 3D topolog- 
ical insulator in classes AIII, Dili, and CI, respectively, 
whereas those in BL classes 2, 7 and 8 can be realized 
either directly on a 2D lattice or, at a surface of topolog- 
ically trivial 3D insulators. Finally we will show in the 
next section by constructing an explicit example, that for 
the BL symmetry class 9_ (CII), the case with minimal 
flavors 7V™ m = 2 can be realized as a surface state of a 
3D topological insulator. 

We will now argue that the gapless nature of 2D Dirac 
Hamiltonians that can be realized at a boundary of a 
topologically non-trivial 3D insulator is stable against 
perturbations V± and A±. 

We first look at the BL classes 1,5, and 6 which corre- 
spond to AIII, Dili, and CI, respectively. These classes 
are special in that the potentials V± are not allowedEJ by 
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the symmetries in these classes, so that 

x=U° + A_- iS o +A+ )- W 

Since the only allowed perturbations are of gauge type, 
one would expect that these perturbations do not spoil 
the gaplessness of the free Dirac spectrum. To see this, 
let us now try to find zero-energy modes. We thus look 
for the solution 

(*_v t+ r+)(£)-> <«) 

where k± = k x ± ik y . We assume that A± have been 
brought into diagonal form by gauge transformations g±, 
so that g±A±g± 1 = A±, where A + = A*_, g\_ = gZ 1 , and 
A± are diagonal matrices with complex entries. Thus, 
the Schrodinger equation for the zero modes reduces to 
0+ + K)XB = {k- + A*+)xa = 0, where X b = g+ Y XB 
and xa = gZ XA- This wave equation has a non-trivial 
solution {xa 7^ or xb 0), only when one of com- 
plex eigenvalues of A + is equal to — k + . The solution 
Xa and xb can then be obtained from xa and xb- It 
should be normalized as x\xa + X B X B = Xa9-9-Xb + 

Xb9+9+Xb = l - ( Notc that 9+ + 9+ 1 )- Thus we con- 
clude that the gaplessness of the spectrum of the 2D 
Dirac modes appearing at the surface of 3D topological 
insulators in AZ symmetry classes AIII, Dili, and CI, 
are robust against arbitrary static perturbations. The 
location of the Dirac cones, however, might be shifted by 
the gauge type perturbations. [As we will discuss, the 
gaplessness of the same Dirac surface modes is also ro- 
bust against Anderson localization arising from random 
perturbations (which break translational invariance), be- 
cause of the presence of Wess-Zumino-Witten (WZW) 
terms.] 

We note that, on the other hand, it is also easy to 
see that for all other classes except All and CII, the 2D 
surface modes have in general a massive spectrum. To 
see this, take V+ = A + = A- = and consider the 
Hamiltonian 

H = k x a x + kyOy + V^a z . (42) 

The potential V- can be diagonalized by a unitary ma- 
trix U, V- -> WV-U =: T = diag(u n ), H -► k x a x + 
kyOy + Ter 2 , and hence the eigenvalues are E n {k x ^ky^ — 

In contrast, for classes All and CII, the potential V- 
satisfies 

V_ = -Vl. (43) 

This guarantees that the corresponding Dirac Hamilto- 
nian has at least one zero eigenvalue when the dimen- 
sionality of V- is odd. We thus cannot completely gap 
out the spectrum by 



V. 3D DIRAC HAMILTONIANS 

The purpose of this section is to construct examples 
of 3D topological insulators in the continuum, one for 
each of the five classes All, Dili, AIII, CI and CII. The 
examples we give are continuum 3D Dirac Hamiltonians 
perturbed by a mass term of some sort. For these exam- 
ples we will compute the winding number v introduced 
in Eq. (|3^) and explicitly derive the surface Dirac modes, 
thereby illustrating the above-mentioned connection be- 
tween the topological properties of the bulk and the exis- 
tence of stable, massless surface states. The Dirac Hamil- 
tonians with the minimal number of components have 
four components for classes All, Dili and AIII, whereas 
the minimal number of components is eight for classes CI 
and CII. 



A. 3D four-component Dirac Hamiltonian 

1. Hamiltonian and its symmetries 

Let us consider the following four-component (3+l)D 
massive Dirac Hamiltonian 

H= —idftOtft + mP, fj, = x,y,z, (44) 

where m G K, and we use the standard (or Dirac) repre- 
sentation of the (3+l)D gamma matrices, 

«M = T»®*M= ^ )> / ? = ^= (J -i)' 

7 5 =Tx=^q^, fi = x,y,z. (45) 

In momentum space, 

H (k) = a,k, + mP = ( k m a ^ ) , (46) 

and the energy spectrum is given by E{k) = 
±V^ 2 + m 2 =: ±A(fc) (two-fold degenerate for each fc). 

At this stage, we have not yet identified the symmetry 
class to which the Dirac Hamiltonian (^) belongs; the 
interpretation of the two gradings represented by a pair of 
standard Pauli matrices, er M and r M needs to be specified. 
As we will see, the four-component Dirac Hamiltonian 
( ff~I| ) realizes topological insulators in classes All, AIII, 
and DHL 

Class All (symplectic) As discussed by Bernevig 
and ChcnjEI the 3D Dirac Hamiltonian (f46| ) is a 
topological insulator in class All since it satisfies 
iUyTL* (k)(—i(7y) = 7i(—k), which we can interpret as a 
TRS for half-integer spin. 

Class Dili In addition to TRS, the Dirac Hamilto- 
nian also satisfies PHS, r y <g> <jyH*{k)T y ®a y = —H(—h). 
Since (r y (g> c y ) T = T y <E> <j yi the 3D Dirac insulator 
( p4| ) can be thought of as a member of class DHL It 
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is possible to unitary transform the Hamiltonian 
H — > diag (<t , —icfy) Wdiag (er , +iu y ) , yielding 



H(k) = 



m 

(-i(Ty)k ■ a 



k ■ a(ia y ) 



M) by 



(47) 



such that the PHS takes on the canonical form displayed 
in Eq. (^). One can easily check that ia y Ti*(k)(~ia y ) = 
H(-k) and T x H(k)r x = -H*(-k). This topological in- 
sulator, Eq. ([47|), describes the fermionic BdG quasipar- 
ticles in the BW state realized in the B phase of liquid 
3 He, for which the d- vector is parallel to the momentum, 
dfe oc k. The /c-dependence of the single particle disper- 
sion Sk of 3 He [see Eg. (|l3|)l is weaker around k = as 
compared to that of the <i-vector and hence neglected. 
I.e., the Dirac mass m here is given by tic minus of the 
chemical potential sf, m — Sk=o = — £f-e3 

Class AIII The 3D Dirac Hamiltonian (Q can be 
viewed as,-an insulator in class AIII due to the (chiral) 
symmetryEa T y TiT y = —Ti. We can bring the 3D Dirac 
Hamiltonian (W4|) into block off-diagonal form by a rota- 
tion Ty — > t z , which transforms the Dirac mass term in 
Eq. (J44|) into the chiral mass term, 



H 



li = x,y,z. 



In momentum space, 
H(k) ■■ 



k-a- 
k • a + im 



(48) 



(49) 



The chiral symmetry is imposed by (3TL(k)f3 — —7i(k). 



2. Wave functions, projector and winding number 

It is well known that the 2D two-component massive 
Dirac Hamiltonian, 'H : ^ TliC (k x , k y ) = k x a x + k y a y + mer 2 , 
is the . sim plest example of a topological insulator in 
2 T&m-. it is an IQH insulator characterized by the 
non-trivial Chern integer (cr xy ), a xy = sgn(m)/2x (e 2 /h). 
If (k x ,k y ,m) is viewed as a set of parameters that can 
be changed adiabatically, the 2D massive Dirac Hamilto- 
nian is nothing but the 2x2 Hamiltonian considered by 
Berry .himself to illustrate the Abelian geometric (Berry) 
phase£J As described below, the 4x4 Dirac Hamiltonian 
( [46| ) can be thought of as a natural generalization of the 
2x2 example Tj 2 ° lac { !<, x h ) aryj ^exhibits a non-trivial 
non-Abelian Berry phase. E§l3l3l3L 



In particular, the two eigenfunctions of the Hamilto- 
nian in Eq. ( |46| ) at wavevector k with negative energy 
E(k) = — A(fc) are given by 



M*0> = 



Mfc)) = 



i 



^2A(A 




(50) 



whereas the eigenfunctions with positive energy E(k) 
+X(k) are 



\Mk)) 



v / 2\(\ - m) 



\Mk)) = 



y / 2X(X - m) 



-k z 


A — m 

k z 
K 
A — m 




(51) 



where k± = k x ± ik y . Note that if m > 0, \u3^(k)} are 
not well-defined at A(fc) = m (i.c, at k = 0). 

Correspondingly, the projector (Q-matrix) onto the 
lowest two negative energy states is given by 



Q(k) = 2P(k) - 1 = — (k^ 



i[3) 



(52) 



By using the chiral grading [see Eq. (|48|)1, we can define 
the g-matrix for class AIII Dirac Hamiltonian (B9h as 



q(k) = -j- (fc M cr M 



(53) 



[Similarly, for the class Dili massive Dirac Hamiltonian 
(p7|), the projector is given by q(k) — ia y (k^a^ — im) / A, 
which satisfies q 1 '(—k) = q(k), in the basis that makes the 
Hamiltonian block off-diagonal as discussed in Eq. ([l2|).] 
The winding number v [Eq. (|3rj)1 for the map represented 
by q(k) from S 3 to U(2) can be computed as 



r i 1 m 
V[q] = 2 ImP 



(54) 



The appearance of a half-integer value for v is common to 
the continuum descriptions, and must be supplemented 
by information about the structure of wave functions at 
high energy (located away from the Dirac point in the 
BZ). See, e.g., the discussion of this issue by Haldane in 
the context of the IQHEE3] 

xupied bands, we can introduce a 



For the lower two-a 



U(2) gauge field byEi 



Af{k)dk, 



( U& (k)\duAk)), a,b 



1,2, 



(55) 



which can be decomposed into U(l) (a ) and SU(2) 
(pP =x > y > z ) parts as 



AJk) = a°(k)— + ai(k)- 



(56) 



While the U(l) part is trivial, the SU(2) part is given by 



a)(k) 



A(A + m) ' 



(57) 



where i = x,y, z and j, I = x,y, z. We have flipped the 
sign of k x , k x — > — k x for notational convenience. 
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3. Boundary Dirac fermions 

We have mentioned above that, quite generally, the 
non-trivial topological properties of the bulk wave func- 
tion manifest themselves as a gapless surface state when 
we terminate the 3D bulk sample by a 2D boundary. To 
see this explicitly, let us take the mass term to be z- 
dependent (m > 0), 



m(z) 



+m, z — > +co, 



-00, 



(58) 



and look for 2D pirae-fjeMBipn solutions localized at the 
boundary z = O.eBMilltidSo For convenience, we take 
the following representation of the massive 3D Dirac 
Hamiltonian in class AIII, H = —ia^d^ — iP r y 5 m(z) . The 
solution to the 3D Dirac equation with energy E(k±) is 



*(*) 




- f z dz m{z ) 



(59) 



where fcj_ = (k x ,k y ) and x± = (x,y) represent the mo- 
mentum and coordinates along the surface, respectively, 
and a(k±) and b(k±) are obtained from the solution to 
the 2D Dirac equation, 



( a{k±) 



V2 



oiarg fc 4 
±1 



(60) 



with E(k±) = ±^jk 2 + ky, respectively. 

Below, we will study the stability of this gapless bound- 
ary 2D Dirac state against the opening of a gap, by per- 
turbing the Hamiltonian by static and homogeneous po- 
tentials which respect the discrete symmetries defining 
the respective symmetry classes. 

Classes All and Dili The gapless nature of the sin- 
gle surface Dirac fermion is protected by TRS since the 
opening of a gap would violate Kramers theorem. In- 
deed, for class All, the only spatially homogeneous per- 
turbation compatible with the TRS is a constant scalar 
potential V, 



+ V, fi = x,y, 



(61) 



which is known not to open a gap. For class Dili, on the 
other hand, even the scalar potential (chemical potential) 
V is prohibited because of PHS. 

The stability of the gapless nature of the single surface 
Dirac fermion is guaranteed by the bulk ^2r4 pYar . ian t in 
the symplectic symmetry class (class All) .EHi3'Eilli3 Al- 
though this protection of the gapless spectrum by the Z2 
invariant also extends to a surface Dirac fermion (which is 
actually Majorana because of PHS in the BdG equation) 
in class Dili when the number of surface Dirac (Majo- 
rana) fermions is odd, it is only the non-trivial winding 
number v in class Dili that guarantees the stability of 



an arbitrary number of gapless surface Dirac (Majorana) 
fermions against perturbations (uniform and random). 

Class AIII A single flavor of 2D Dirac fermions in 
class AIII can be perturbed by a static and homogeneous 
vector potential: 



H 



11, V = x,y. 



(62) 



The vector potential perturbation shifts the location of 
the node, but does not open a gap. 

Although there is a "hidden" TRS in class AIII, the 
stability of this single Dirac fermion ( |62| ) is not protected 
by the Z2 invariant. This is so since in order to reveal the 
TRS, we need to consider the full BdG Hamiltonian H4 in 
(||); rather than H2 defined in ( |l6| ) and (p6[). In the 
number of flavors of the surface Dirac fermions is counted 
as two, and not protected by the Z2 invariant. Again, it 
is the winding number v that guarantees the stability of 
an arbitrary number of flavors of gapless surface Dirac 
fermions against perturbations. 



B. 3D 8-component Dirac Hamiltonian 

It turns out that in general we cannot have a (3+l)D 
4-component Dirac Hamiltonian which is a member of 
class CI/CII, and which also possesses a gapless Dirac 
fermion surface mode. We are thus led to consider a 
(3+l)D 8-component Dirac Hamiltonian. (It is possible 
to construct gapless four-component Dirac Hamiltonian 
in classes CI and CII, but we cannot give a mass to them.) 

Class CI The following massive 3D 8-componcnt 
Dirac Hamiltonian 



n 



o 



D 




D = io- y (3 



17717 



(63) 



is a member of class CI since D T (k) = D(—k). [See Eq. 
(p5|).] The energy spectrum at wavevector k is given by 
E{k) = ±\/k 2 + m 2 = ±A(fc), where each eigenvalue is 
four-fold degenerate. 

The projector takes on block off-diagonal form and is 
given by 



Q(k)=2P(k)-l = ~H(k), 



q{k) 



-^ivyP (k^a^ ~ irwy 5 ) 



The winding number can be computed as 



r n 1 m 

2 771 



(64) 



(65) 



which is twice as large as the winding number for the 
four-component case. As before, this winding number 
should be interpreted cither m/|m| x 2 or 0, depending 
on the behavior of the wave function at higher energy. 
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When we terminate the 3D Dirac insulator (p3|), by a 
2D boundary, by making the mass term z-dependent as 
in Eq. (pg), we find two flavors of surface Dirac fermions, 



n 

D 



D 
£>t 

Ty(k^ 



A. r a x 



A y a y + A z a z ) , 



(66) 
Z al- 



wherc we have included perturbations A x ^ y _ z € 
lowed by class CI symmetries. One can easily check 
that D T (k) = D(—k). The gaplcss nature of this four- 
component Dirac fermion is stable against arbitrary val- 
ues of 3 complex (6 real) parameters A Xj y jZ . Indeed, 
just like a vector potential perturbation in class AIII, 
shifts the location of the Dirac node from (0, 0) to 



(k x , ky), where (k x , fc°) is a solution to 
(k° x ) 2 - (k° v ) 2 - RcA 2 +lmA 2 



0, 



k x k° + (ReA ■ ImA) = 0. 



(67) 



Class CII The following 8-componcnt 3D Dirac 
Hamiltonian, 



H = 



D 
£>t 



D = k. 



+ m(3 = D\ (68) 



is a member of class CII since ia y D*(k)(—ia y ) = D(—k), 
and has a gapped spectrum. E(k) — ±\/k 2 + m 2 . The 
projector Q and g-matrix are given by 



Q(fc) = 2P(fc) - 1 = -\n{k), 

q(k) = -- (k^a^ + m/3) . 



(69) 



Observe that, compared with the 3D Dirac insulator 
in class CI, Eq. (|63|), the mass term for the 3D Dirac 
insulator in class CII, Eq. (|68|), is given by the Dirac 
mass term (m/3), not by the chiral mass term (irwy 5 ). 
Due to this difference, the winding number vanished for 
the 3D Dirac Hamiltonian in class CII, Eq. (|6li|): 



'[?] = o. 



(70) 



In spite of the vanishing of the winding number, we do 
find two flavors of two-component Dirac fermions at the 
surface of a 3D Dirac insulator in class CII (|68), when 
making the mass term z-dependent as in Eq^5q). In 
particular, consider the genera l form of the Dirac Hamil- 
tonian on the 2D surfaccE3li^ 



n 

D 



D 

Dt / 

v + fc_ + a 
fc + — a v*. 



(71) 



= (k x + ia x ) o~ x + (k y + ia y ) a y 
+ Reu+do + Hmv + (T Zl 

where the perturbations a x . y and v + are the only ones 
allowed by class CII symmetries. (We used the notation 



a = Re a + ilma = a y + ia x .) It turns out that the gap- 
less nature of these surface Dirac fermions is preserved 
by these perturbations. To see this, consider the deter- 
minant of the Hamiltonian 



det (£>£>*) = (K| 2 - \k+\ 2 + \a\ 2 f - (afc_ - ak+) 2 , 



(72) 



which vanishes when 



l fc +| 2 = K'+l 2 + M 2 and (k x ,k y ) _L (a x ,a y ). (73) 

This shows that it is always possible to find a wavevector 
(k x ,k y ) for which the determinant, and thus the energy 
eigenvalue vanishes, proving the absence of a gap. 

Therefore we conclude that the 3D Dirac insulator ( |6"8| ) 
is a non-trivial topological insulator in class CII: it is im- 
possible to deform the insulator ( |S8[ ) into a topologically 
trivial insulator (an insulator without a stable surface 
state) without closing the energy gap in the 3D bulk, be- 
cause the existence of the gapless Dirac fermion surface 
modes plays the role of a topological invariant. On the 
other hand, when the number of flavors is twice an even 
integer, one can easily find a perturbation that gives a 
mass gap to all surface Dirac fermions. 



VI. TOPOLOGICAL FIELD THEORY 
DESCRIPTION 

In order to understand more intuitively the reason why 
symmetry classes with a sublattice (chiral) symmetry 
(classes AIII, Dili, and CI) possess stable gapless surface 
Dirac fermion modes, we derive in this section a doubled- 
Chcrn-Simons field theory describing the 3D bulk insula- 
tor. To this-end we identify, following the spirit of Read 
and GreenJ^I conserved charges of the action, and intro- 
duce external gauge fields that couple minimally to these 
charges. The gapped fcrmionic degrees of freedom in the 
3D bulk are then integrated out to derive the effective 
action of the gauge fields. 

A similar procedure has also been discussed for 
domain-wall fcrmionsLj in lattice gauge theory, where a 
(non-doubled) Chern-Simons theory can be derived lot 
[(2n- 1)+1]D boundary fermions of the (2n+ 1 )D bulk.EE 

As an example, let us take the class AIII Dirac insu- 
lator ( [l8| ) in three spatial dimensions. The generating 
function for the single-particle Green's function can be 
written as a fermionic functional integral, 



S = i d 3 x^ (H-irj)il). 



(74) 



Here, note that we are using a three-dimensional (Eu- 
clidean) action, instead of a (3+l)-dimcnsional one (com- 
pare, e.g., with Rcf. 16), since we are focusing on single- 



particle Green's functions in the absence of interactions. 



16 



A finite lcvcl-broadcning term ?/ ^ is necessary to reg- 
ularize delta functions appearing in the single particle 
Green's function. 

The action enjoys a (electromagnetic) U(l) symmetry 



tp' — > tjj'e 



1> 



-il) 



(75) 



where 9 6 M. Due to the sublatticc (or chiral) symme- 
try, /3W/3 = ~TC, the functional integral possesses the 
additional symmetry 



0, 



(76) 



when 77 = 0. A finite level-broadening term 77 spoils this 
symmetry. Instead of adding a level-broadening term, 
however, regularization of the functional integral can also 
be achieved, alternatively, by attaching ideal leads (or 
perfect absorbers) respecting the chiral symmetry to the 
sample. 

We now proceed to derive the Chern-Simons theory. 
Corresponding to the two continuous symmetries dis- 
cussed above, we couple two U(l) gauge fields a M and 



to the fermions, 



d 3 x ip(@-ifi- ij $ + 77175) tp, (77) 



where we have introduced the abbreviations xp = ip'P, 
7° = /?, 7 fe = (3ak, and $ '■= J^a-n- While the ex- 
ternal U(l) gauge field Ou, associated with the global 
U(l) symmetry of Eq. (f75|), couples to the electromag- 
netic current, the external "axial" gauge field detects 
the sublattice-resolved current, as it is associated with 
the global U(l) symmetry transformation defined in Eq. 
(f76|), where equal and opposite U(l) transformations are 
performed on the two sublattices A and B of the under- 
lying bipartite lattice. 

To be more general, we discuss the case of TV repli- 
cas of the above 3D Dirac fermions, and couple them 
minimally with two U(2V) gauge fields, a M = a^T a and 
fop = b^T a , with the generators T a . The use of repli- 
cas is a convenient method to compute disorder averaged 
physical quantities in the presence of random impurities. 
Since we only intend to give a schematic derivation of the 
doubled Chern-Simons theory, we do not add any explicit 
disorder potential. 

We now integrate out the fermions and derive the ef- 
fective action for the gauge fields a M and 6 M , 



-Soft [a^,6 M ] 



(78) 



by a derivative expansion 

S eS = -Trln(Go 1 -V) 

OO 1 

= -Tr lnG^ 1 + V -Tr {G V) n , 



(79) 



where Go denotes the propagator of free 3D Dirac 
fermions, which is given in momentum space by 



Go(fe) 



i# + 77175 



777/ 



whilst 

V{q)^-i^{q)-i 1 ^{q). 
Introducing the linear combinations 



(80) 

(81) 
(82) 



the resultant effective action, to leading order in the 
derivative expansion, takes the form of a (Euclidean) 
doubled Chern-Simons theory, 



1 777 

S eS = -{?-(l[A+)-I[A-])+dW. 

2 m v ' 



(83) 



2/ 



where "div." represents an ultraviolet (UV) linearly di- 
vergent piece. This divergence is closely related to the 
appearance of the half-integer coefficient of the Chern- 
Simons term, sgn(m) x 1/2: the action ±sgn(ra)/L4]/2 
is not gauge invariant by itself. [See also the discussion 
below Eq. (U).] 

Introducing a gauge invariant regulator, such as the 
Pauli-Villars (PV) regularization, cures, both .the UV di- 
vergence and the half-integer coefficient.E23E23 Here, note 
that, since there is no chiral anomaly in 3D, the func- 
tional integral can be regularized without breaking the 
two U(iV) gauge symmetries, although the parity symme- 
try can be destroyed by the regularization. In the Pauli- 
Villars regularization, we define the physical, divergence- 
free effective action S^g by 



c-PV 
D cff 



S off (0)- lim S off (M), 

M 2 — *oo 



(84) 



where S e g(M) represents the effective action in the pres- 
ence of two massive Dirac particles: the original particle 
with the mass m and another one with mass M, which 
we take to infinity (M 2 — > 00). Here, the second par- 
ticle (which is bosonic) might be interpreted as supple- 
menting the missing information far away from the Dirac 
point discussed around Eq. (|54|). The coefficient of the 
Chern-Simons terms /L4 + ] and I[.A _ ] in , which is 
(1/2)(tt7/H - M/\M\) instead of (l/2)m/|ra|, depends 
on the sign of the regulator mass M: when sgn(M) = 
— sgn(m), this coefficient equals sgn(TO) whereas it van- 
ishes when sgn(M) = +sgn(?77). These two cases rep- 
resent the topological non-trivial and trivial phases, re- 
spectively. 

Once we have established the appearance of the dou- 
bled Chern-Simons term for the resulting 3D bulk theory, 
we conclude that the surface degrees of freedom, which 
appear when the 3D bulk is terminated by a 2D sur- 
face, are described by the two-dimensional U(iV) Wcss- 
Zumino-Witten (WZW) theory at level k = 1. The 
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WZW theory is well known to be gapless and tp-oossess 
both, holomorphic and antiholomorphic sectors.E23 Thus, 
the U(iV) x V(N) symmetry of the gapped 3D bulk the- 
ory, which represents two independent transformations 
for each sublattice, turns into the two independent holo- 
morphic and antiholomorphic U(-ZV) gauge symmetries of 
the WZW theory describing the resulting degrees of free- 
dom at the surface. 

An entirely analogous discussion can be carried 
through for the 3D topological insulators in symmetry 
classes Dili and CI, for which the relevant gauge groups 
appearing in the gapless WZW theory at the 2D surface 
are 0(2iV) and Sp(iV), respectively. 



VII. MANY-BODY ASPECTS OF 3D 
TOPOLOGICAL SUPERCONDUCTORS 

Up to now, we have treated the three-dimensional su- 
perconductors at the mean-field level of pairing (where 
the gap function A is "frozen" to be a constant, i.e., 
frequency independent): we have focused solely on the 
dynamics of the fermionic BdG quasiparticles, existing 
deep within the superconducting state. In the full de- 
scription of superconductivity however, the gap function 
A, together with the electromagnetic U(l) gauge field, 
have to be regarded as dynamical entities. Once we em- 
ploy this full description, the superconductors of interest 
in this paper should be regarded as topological phases 
with non-trivial pmind state degeneracy (and charge 
fractionalization) j£3l£j (Specifically, while the BdG quasi- 
particle sector of the B-phase 3 He is a topological Dirac 
insulator in class Dili [as discussed in Eq. (f47|)1, the B- 
phase of liquid 3 He is in reality not a 3D topological 
phase, as it supports gapless Nambu-Goldstonc bosons.) 



A. Many-body wave function in real space 

Let us first take a closer look at the many-body 
ground state wave function of the three-dimensional 
class Dili Dirac insulator [see e.g., Eq. (E^)]. In a 
triplet superconductor, the BCS ground state |BCS) 
projected onto a space of fixed electron number N, 
\x\, <J\]X2, &2] ■ ■ ■ ; xn, on) (where Oi represents spin co- 
ordinate), is given by the wave function 

*(.ti, £7i ; x 2 , cr 2 ;...; x N , <J N ) 

■ = (xi,ai;x 2 ,cr 2 ; . . . ; x N ,a N \BCS) 

= Pf \ff(xi,(Ti;Xj,(Tj)] , (85) 

where Pf denotes the Pfaffian of the matrix gij := 
g(xi, <7i\ Xj, cTj ) .til The Fourier transform of g(x, a, y, r) = 
g ar {x — y) as obtained from j5C| ) and (|5l|), reads 



,{k) = {-X + m) 2fc2 » . 



Noting that \(k) = ^k 2 + m 2 -» |m| + |fc| 2 /(2|m|) H 

in the long- wavelength limit, k — ► 0, the expression in 
Eq. (pq) takes in that limit the form 



9(k) 



(k ■ a)ia y 
4|to| ' 

(k- a)ia y 
H — ^ JL - 



m > 0, 



m < 0. 



(87) 



Correspondingly, the real-space wave function g(r) = 
(27r) -3 / d 3 k e zk ' r g(k) takes at long scales the following 
form 



a a ia. 



id^ 3) (r), m>0, 



i\m\ " 
11 47rr 3 



m < 0. 



This behavior is similar to the strong and weak pair- 
ing phases of JLhc (2+l)-dimensional chiral p-w&ve 
superconductor .Eil In one phase, the strong pairing phase 
(to > 0), the wave function g(r) of a pair is short-ranged, 
whereas in the other, the weak pairing phase (to < 0), 
g(r) exhibits a power-law behavior and is given by the 
correlation function of the two-component massless 3D 
Dirac fermion. Thus, in the weak pairing phase, the 
many-body ground state wave function behaves at large 
scales as 



W(xi, £7i ; x 2 , cr 2 ; . . . ; x N , cr N ) 

'(cr- (xi-x^iay)^ 



(89) 



Observe that this is nothing but the multi-point corre- 
lation function of a (simple) 3D conformal field theory, 
namely the 3D free Majorana fermion quantum field the- 
ory defined by the partition function 



■ Jd 3 xC 



(90) 



(86) 



where ip is a two-component Grassmann variable with 
Majorana condition -ip — ip T i<r y , and /i = x,y,z. This 
is analogous to the Moore-Read Pfaffian wave function, 
which is given by the multi-point correlation function 
of the 2D Ising conformal field theory (free Majorana 
fermion field theory). 



B. Ground state degeneracy 

With both the pairing potential A and the U(l) gauge 
field being frozen, there is a unique ground state both 
in the strong and weak pairing phases. We now include 
quantum fluctuations of A and the U(l) gauge field. One 
consequence of the inclusion of these as dynamical de- 
grees of freedom is the appearance of a non-trivial ground 
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state degeneracy. The counting of ground states in each 
phase is completely parallel to the case of the Moore- 
Read Pfaffian state as we will see below. 

To count the ground state degeneracy on the three- 
torus T 3 , we consider periodic or anti-periodic bound- 
ary conditions (BCs) along the three cycles in the 
x,y,z directions. We denote sectors with different BCs 
by (i x , L y , l z ), where = ± represents periodic/anti- 
periodk, BC. Following the argument by Read and 
Green,Eil we notice that the state with k = is allowed 
only for the (+, +, +) sector, and it is occupied in the 
weak pairing phase whereas it is unoccupied in the strong 
pairing phase. 

In the strong pairing phase, different boundary con- 
ditions lead to 2 3 degenerate ground states, and all of 
them have an even number of fcrmions. On the other 
hand, in the weak pairing phase, the ground state for the 
boundary condition sector (+, +, +) has an odd num- 
ber of particles because of an additional occupied state 
at k = 0. Thus, the ground state degeneracy for an 
even number of fermions is 2 3 — 1 — 7 whereas there is a 
unique ground state for an odd number of fcrmions. This 
should be contrasted with the ground state degeneracy 
of 2 3 present in the 3D Abelian Higgs model which can 
be described by a (3 + 1)D BF topological field theoryS 
It is unclear what kind of bulk topological field theory 
can describe the weak pairing phase, as it has fermionic 
excitations at boundaries, unlike the bosonic boundary 
excitations in the strong pairing phase described by the 
BF topological field theory. 

The smaller topological degeneracy in the weak pairing 
phase can also be understood in terms of the "blocking 
mechanism" introduced in Ref. ^35|. The smaller topologi- 
cal degeneracy in the Moore- Read Pfaffian state happens 
because a vortex-antivortex excitation carries a Majo- 
rana fermion at the core. In a topological phase, a differ- 
ent ground state, starting from a given ground state, can 
be generated by first creating a particle- antiparticle pair 
out of the ground state, then moving around the quasi- 
particle along a homotopically non-trivial cycle, and fi- 
nally pair-annihilating the pairs. If the quasi-particle ac- 
commodates a Majorana fermion, however, the last step 
of the above procedure, which is pair-annihilation, might 
not be possible (it might be "blocked"). In the 3D Pfaf- 
fian state, we have vortex lines, insteadjaf. ^rtires., which 
do support Majorana fermion modes .ooE£ZrE23 Thus, 
we expect a similar blocking mechanism should apply. 

We now briefly discuss the effects of intcrparticlc in- 
teractions. Since short-range interactions are irrelevant 
by power-counting for free Dirac fermions in (2 + 1) di- 
mensions, wc would expect the gaplcss fermionic surface 
modes in the weak pairing phase to be stable against 
the formation of a gap, up to, possibly, some critical in- 
teraction strength (certainly when random disorder po- 
tentials are not simultaneously present). This should be 
contrasted with the surface states in the strong pairing 
phase, which are generically gapped as wc can_sp£j_fjpr ex- 
ample, from the BF topological field theory.lEBL-3 This 



should also be compared with the surface states of three- 
dimensional Z 2 topological insulators in the symplectic 
class (All), which are unstable against the BCS pairing 
instability because there is a finite Fermi surface (circle), 
i.e., finite density of states, within the surface Brillouin 
zone, for a general value of the chemical potential. 



VIII. DISCUSSION 

In this paper, we have undertaken the program of clas- 
sifying possible phases of topological insulators and su- 
perconductors in three spatial dimensions. Our results 
have their root in the very general classification scheme 
for random matrices obtained by Zirnbaucr, and Altland 
and Zirnbauer (AZ)! 3 ^! 3 ^ more than a decade ago, result- 
ing in ten such classes which extend the well-known three 
Wigner-Dyson classes. Guided by the lessons learned 
from [ -jiU|'_Z2 topologicafjasulator discussed by Kane and 
McIcQOEj and othersJalilJ which belongs to the symplec- 
tic ( "spin-orbit" \ symmetry class (All) in the AZ classifi- 
cation scheme,L3E3 we have asked if two quantum ground 
states in a given symmetry class can be continuously de- 
formed into each other while keeping the discrete symme- 
tries defining the symmetry class intact. Specifically, we 
have shown that, in addition to the three-dimensional Z2 
topological insulators there exist 3D topological insula- 
tors possessing the symmetries of four additional random 
matrix classes denoted by AIII, Dili, CI and CII in the 
work of AZ, all of which support stable gaplcss Dirac 
fermion surface modes (Majorana fermion surface modes 
for class Dili). In particular, we find that the topological 
properties of the bulk wave functions in the three sym- 
metry classes AIII, Dili and CI are characterized by an 
integral winding number, while the bulk characteristics of 
topological insulators in class CII can be described by a 
Z2 number, akin to the well studied topological insulator 
in the symplectic symmetry class AIL 



A. Topological bulk characteristics and Anderson 
derealization at the boundary 

Another lesson learned from the Z2 topological insu- 
lator is an intimate connection between the topologi- 
cal characteristics of the clean (no disorder, or impu- 
rities) system in the 3D bulk and the Anderson lo- 
calization physics occurring, due to disorder, at two- 
dimensional boundaries of such a system: the surface 
of a three-dimensional Z2 topological insulator is a per- 
fect metal (Z2 topological metal) in the presence of dis- 
order which respects the TRS. This can be understood 
in terms of the, fifild , theoretical framework of Ander- 
son localization 



loqjiiq 



The fermionic replica non-linear 
a model (NLctM) describing quantum transport in the 
corresponding symplectic ("Wigner-Dyson") symmetry 
class possesses the coset space 0(4iV)/0(2iV) x 0(2N) as 
target space. (N is the number of replicas.) Because the 
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AZ class 


Spa.ce of transfer matrices 


2D top. insulator 


Possible phys. realization 


A 

i\ 


U(J>, q)/U(p) x U(q) 


lytlJi (p 7^ q) 


Vjrd,.rlS / .rllVjTcLrVb 


AT 








All (even) 


SO*f4An/UY2An 






All (odd) 


SO*(47V + 2)/U(27V + l) 


Z 2 top. ins. (QSHE) 


HgTc/(Hg,Cd)Te 


A TTT 








RnT 
















D 


SO„(p,g)/SO(p) x SO(g) 


Thermal QHE (p =fc q) 


Spinless chiral p-wave SC 


C 


Sp(p, 9 )/Sp(p) xSp(g) 


Spin QHE (p / g) 


(ci ± id)-wave SC 


Dili (even) 


SO(2A,C)/SO(2A) 






Dili (odd) 


SO(2iV+l,C)/SO(2JV + l) 


Z 2 top. SC 


(p + ip) x (p — ip)-wave SC 


CI 


Sp(JV,C)/Sp(JV) 







TABLE IV: This table lists the space of ensembles of (the radial ee«rdinates of) transfer matrices for quasi one-dimensional 
disordered quantum wires for each Altland-Zirnbauer (AZ) class.lii3 Five of these ensembles of transfer matrices describe 
localization properties of an edge of a two-dimensional topological insulator or superconductor (SC). The conventional name 
of these five two-dimensional topological insulators (2D top. ins.) is given in the third column. The last column lists some 
possible physical realizations of these topological insulators. 



homotopy group 7r 2 [0(4A^)/0(2iV) x 0(2iV)] =Z 2 ,aZ 2 
topological term is allowed in the action of this NLerM, 
and this term is indeed realized at a surface of 3D Z 2 
topological insulators. It is the Z 2 topological term that 
is responsible for the lack-pih^alization and the metallic 
behavior at the surface .BEjE§BEl 

A key result of the present paper is a generaliza- 
tion of these properties of the 3D Z 2 topological insu- 
lator to topological insulators belonging to the above 
mentioned four symmetry classes AIII, Dili, CI and 
CII. For three of these four classes, namely for AIII, 
Dili, and CI, which describe the dynamics of quasipar- 
ticles within certain superconductors, the corresponding 
NLerM describing Anderson localization at the surface of 
the 3D bulk is the principal chiral model (PCM) on the 
groups V(N), 0(2N), and Sp(A^), respectively, supple- 
mented by a Wess-Zumino-Witten (WZW) term. For 
class CII, the corresponding NLoM is denned on the 
coset space XJ(2N)/0(2N), which allows for a Z 2 topo- 
logical term, since the homotop y g roup of this space is 
7r 2 [U(2iV)/O(2A0] = z 2- Table M summarizes the tar- 



get spaces of the corresponding fermionic replica NLaMs, 
as well as possible 2D topological or WZW terms. (See, 
e.g., Refs. ^ and ^l]). The NLctMs living on the surface 
of the 3D bulk "remember" the non-trivial topological 
characteristics of the bulk througk-tbe presence of either 
a WZW or a Z 2 topological termJlB Due to these WZW 
or topological terms, the quantum states existing at the 
surface of the 3D topological bulk arc gaplcss and their 
gaplessncss is protected against (Anderson) localization 
by random potentials respecting the discrete symmetries. 
We thus conclude that surfaces of 3D topological insula- 
tors with a TRS are always ( "topologically" ) delocalized. 
[For sublattice and superconducting classes, the corre- 
sponding gapless states at the surface are semi-metal like 
(with conductivities of order unity in natural units; see 
end of this section), unlike the surface states of Z 2 topo- 



logical insulators, which are perfect metals. 



Topological insulators in one and two spatial 
dimensions 



This correspondence between nontrivial topological 
characteristics of an insulator in the bulk and lack of 
Anderson localization due to random impurities at the 
boundary also applies to 2D topological insulators and 
their ID edge modes at boundaries. In (quasi) ID, An- 
derson localization problems can be well-described by 
the Dorokhov-Mcllo-Pcreyra-Kumar (DMPK) equations, 
which is a Fokkcr-Planck equation describing the distri- 
bution of the eigenvalues of transfer matrices of the quasi 
ID wire as a function of the wire length. The ensembles of 
transfer matrices can be systematically enum erated , and 
there are twelve possible DMPK equations.! 112 !! 113 ! (See 
Table [V.) The extra two cases which appear here, in 



addition to the ten AZ symmetry classes, cannot be re- 
alized as a quasi ID tight-binding (lattice) model, but 
can only be realized at the (one-dimensional) boundary 
of 2D topological insulators. One ofpthese two cases is 
the symplectic symmetry class (AII),EjEj which can be 
realized, for example, at the edge o£i;he 2D Kane-Mele 
model. The other is in class DHI,E-3 and can be real- 
ized, for example, at the edge of the equal superposi- 
tion of two chiral p-wave superconductors with opposite 
chiralitics [faj+.jm)-. and (p — ip)-wave], in two spatial 
dimensions.HOOlI 2 ! In classes A (unitary), D, and C, the 
DMPK equation depends on two integers p and q, rep- 
resenting the number of left- and right-moving channels, 
respectively. When the numbers of left- and right-moving 
channels are not equal (p ^ q), i.e., when the quasi ID 
system is chiral, the corresponding Anderson localization 
problem can only be realized at an edge of a 2D topolog- 
ical insulator. Specifically, in class A this topologically 
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AZ class 


Space of Hamiltonians 


ID top. ins. 


A 






AI 


IK AO/Of AO 




All 


\](2N)/Sp(N) 




A TTT 


v(P + q)/ v(p) x U(q) 


iLi 




o^Jyp -+- q) / d*U\p) X o\J\q) 


7/ 


PIT 


op(p + q)/z>p(p) x ap^gj 


1> 
in 


D (even) 


SO(N) x SO(Ar)/SO(./V) 


- 


D (odd) ("B") 


SO(2A + 1) 


Z 2 


C 


Sp(JV) x Sp(JV)/Sp(AT) 




Dili (even) 


SO(2iV)/U(JV) 




Dili (odd) 


SO(4A + 2)/U(2A + 1) 


z 2 


CI 


Sp(iV)/U(iV) 





TABLE V: This table lists the space of ensembles of ran- 
dom matrices for zero-dimensional disorder 



quantum sys- 
tems for each Altland-Zirnbauer (AZ) class.EiiS According to 
Ref. 115 , five of these ensembles of random matrices describe 

(zero-dimensional) boundary of a 



localization properties at 
one- dimensional topological insulator (ID top. ins 



non-trivial 2D quantum ground state is commonly known 
as the IQHE, in class D it_is_thc thermal quantum Hall 
effect in superconductors,EJu and in . clas s C the spin 
quantum Hall effect in superconductors JZjO respectively. 

Finally, a similar correspondence exists also between 
one-dimensional topological insulators and their zero- 
dimensional edges, where disorder effects can be dis- 
cussed in terros-pf random matrix theories (RMTs). Fol- 
lowing IvanovElj, five out of ten AZ classes, Dili, D, AIII, 
BDI, and CII, allow random matrix ensembles with ex- 
act zero eigenvalue (s), or "zero- modes". (See Table 0.) 
This in turn suggests the existence .aLonc-dimcnsional 
topological insulators in these classesti£l. To summarize, 
we have listed all topological insulators as a function of 
symmetry and spatial dimension in Table |. The labels Z 
and I12 used in this table indicate whether the different 
topological sectors can be labeled by an integer or a Z2 
quantity, respectively. 



C. Experimental implications 

The topological insulators discussed in this paper can 
be realized in nature: as explained in Sec. VII, the 
fermionic sector of the quasiparticles in the B-phase of 
liquid 3 He is an example of the topological insulator in 
the superconductor class DHL Also, an unconventional 
superconductor in a heavy-fermion compound, say, could 
be a possible realization of a 3D topological phase, pos- 
sessing, e.g., a non-trivial ground state degeneracy. An- 
other arena where these novel 3D quantum states with 
exotic topological characteristics may be realized exper- 
imentally is that of cold atom systems (tunable via the 
p-wave Feshbach resonance). [The realization of a 2D 
topological phase (the Moore-Read Pfaffian state or chi- 



ral 2D p-wavc superconductor) in co ld a tom systems has 
recently been discussed in Rcfs. 117 and 118.] Unlike the 
B phase of liquid 3 He, such a realization might allow us to 
go back and forth between the weak and strong pairing 
phases of the topological insulators by detuning. Last 
but not the least, strong correlations among electrons 
(or spins) might spontaneously give rise to these topolog- 
ical phases,tij by forming a non-trivial band structure 



for fermionic excitations (e.g., spinons), which can be ex- 
plorcd-bjf, e.g., slave-particle mean-field theories of a spin 
liquidJl23 

One of the direct signatures of the non-trivial topo- 
logical characteristics of the quantum state in the three- 
dimensional bulk is the appearance of gapless relativistic 
fcrmion modes at surfaces terminating the bulk, which 
are stable to interactions and to disorder. It should 
be possible to detect these surface states using vari- 
ous experimental probes, such as tunncling/STM probes, 
and especially angle resolved photocmission spectroscopy 
as already done successfully in the Bismuth- Antimony 
alloys Ea 

Of particular interest are transport measurements, as 
the 3D topological insulators always possess delocalizcd 
gapless modes propagating at their surface, even in the 
presence of disorder. For the 3D Z2 topological insulators 
such surface modes (in the symplertic .symmetry class) 
are predicted to be a perfect metal.OEZI Hence electrical 
transport measurements can be used to determine if a 
specific insulating material is a Z2 topological insulator 
or not, which is a test independent of, say, photocmission 
experiments, in which one counts the number of surface 
Dirac fermion flavors. 

We suggest that, for the superconducting classes DHL 
AIII, and CI, it would be interesting to measure cither 
spin transport (for spin conserving symmetry classes, 
AIII and CI) or thermal transport (for all the three 
classes) properties of the gapless delocalized surface 
modes. The spin conductivity {<J s xx ln ) or the thermal 
conductivity divided by temperature (k xx /T) is, in the 
absence of disorder, of order unity (in natural units) be- 
cause of the vanishing density of states at zero energy: 
when the Dirac cone is isotropic (i.e., the fermi velocity 
is the same in all directions}, the spin conductivity in 
classes AIII and CI is givena by a x p x ln = 1/tt X s 2 /h per 
Dirac fermion, where s = 1/2 is the spin "charge" car- 
ried by quasiparticles. The thermal conductivity is then 
given by k xx /T = 47r 2 /(3f7^.P in ) valid for all three classes 
AIII, Dili, and CI (where for classes AIII, and CI this 
represents the Wiedemann- Franz law). For each surface 
Majorana fermion in class DHL the thermal conductivity 
is half the value obtained for a single Dirac fermion, and 
the contributions from several flavors are additive. 

The values for these transport coefficients are com- 
pletely robust against disorder. This can be directly ob- 
served for the case of minimal number of surface Dirac 
fermions; Eq. @ with V = for class Dili, Eq. @ for 
class AIII, and Eq. ([36]) for class CI. For class Dill, there 
is simply no disorder potential allowed by symmetries. 



21 



While, as discussed, gauge type randomness is possible in 
class AIII [random U(l) gauge field] and class CI [random 
SU(2) gauge field], itis jsMfpwfo tl^at they do not affect the 



transport coefficient 



Thus, at the surface of 



the 3D topological insulators in class Dili, AIII, and CI, 
the transport coefficients (spin and thermal conductivi- 
ties) are temperature independent and universal. 
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